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Abstract. We have used the optical lines of N II and Fe III to study the wind of the luminous blue variable 
P Cyg. This was performed by applying a version of the Self Absorption Curve (SAC) method, involving few 
assumptions, to lines whose flux can be measured. A rather surprising result was obtained; the lines of more 
excited multiplets without blue shifted absorption components appear to be optically thick, while the lines of the 
most excited multiplets may show some indications of being optically thicker than the lines of less excited ones. 
Explanations of such effects are discussed, including possible inhomogeneities in the wind. 
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1. Introduction 

The observational history of the hypergiant P Cygni be- 
gun about 400 years ago, when it was discovered by 
Willcm Blacu as a nova-like object. Because of its spec- 
tral and photometric characteristics the star is considered 
to be the prototyp e of th e P Cygni- type stars (PCT), as 
defined by Lamers ( 1986 ), and as one of the prototypes of 
the giant eruption Luminous Blue Variables (Humphreys 
1999] ). In addit ion, r ecent studies (Markova et al. 2001 



de Groot et al. |200l| ) showed that over the last 15 years 
PCygni has behaved as a typical weak S Doradus variable 
in a short S Dor phase (van Genderen |200l| ) . More infor- 
mation about this enigmatic object can be found elsewhere 
(Israelian and de Groot |1999|, ASP Conf Ser. 233). 



Stahl et al. (1993) reported the presence of a num- 
ber of pure emission lines in addition to the dominant 
P Cygni-type lines in the optical spectrum of P Cygni. 
Later Markova & de Groot (1997) showed, comparing line- 
identification lists from various observational epochs, that 
the optical emission spectrum of P Cygni was much richer 
and intense in the mid nineties of the nineteenth century, 
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than 60 years ago. In particular, the authors noted that 
more than 70% of the pure emission lines appear to be, 
of recent origin. Among these are forbidden lines of Fe II, 
Ni II, Ti II Fe III and N II; high excitation lines of Fe III 
and N II and lines of low and medium excitation of Si II. 
If not due to observational selection caused by the con- 
tinuous improvement of the signal-to-noise ratio of photo- 
graphic spectra during the last century, this result might 
indicate the presence of a very long-term variation in the 
wind+photospheric properties of the star. In this context, 
it seems worthwhile to try to obtain additional informa- 
tion on the nature and the origin of P Cygni's emission 
spectrum using methods not used up to now, such as for 
instance the SAC (Self- Absorption Curve) method devel- 
oped by Friedjung and Muratorio ( |1987 ). This method in- 
volves a semi-empirical analysis of emission line spectra of 
complex atoms and ions, without assuming detailed mod- 
els for the objects emitting these spectra. By means of the 
SAC method valuable information (e.g. the self absorption 
effects, the level population laws and population anoma- 
lies) has been obtained for the Fe II emission spectrum of 
many stars of differ ent types, such as AG Car (Muratorio 
and Friedjung |"l988|), the VV Cep star KQ Pup (Muratorio 
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et al 1992| ), while procedures for using the method is de- 
scribed by Baratta et al (1998). Use of the SAC method 
can require however a few hypotheses about the relative 
populations of levels inside the same spectroscopic term 
and the distribution of self-absorption relative to that of 
the emission of different lines in different parts of the line 
emitting region. The second of these hypotheses was in 
any case not needed in the present work. 

The main purpose of our study is to apply the SAC 
method to the Fe III and N II emission spectra of P Cygni. 
In this way we can hope to gain additional information 
about the nature and the origin of these spectra and thus 
to get a deeper insight into the physics of the star it- 
self. The observational material is described in Section 2, 
where a number of problems concerning the derivation of 
the net emission equivalent widths of the studied lines are 
also discussed and resolved. Section 3 presents the results 
obtained through the SAC analysis while Section 4 deals 
with the interpretation of these results. 



2. Observations and data reduction. EW 
measurements 

The present analysis is based on one spectrum taken at 
the Heidelberg Observatory on August 9 1991 with a fibre 
echelle spectrograph attached to a telescope of the 70 cm 
class (Stahl et al p995|). An EEV CCD with 1252 X 770 



pixels of 22 /i size was used as a detector. The spectrum 
covers an interval of 2 700A, from A4001 to A6773, with a 
spectral resoluti on of 12 000 and a signal to noise ratio of 
333 (Stahl et al. 1993 ). In cases where measurements were 
doubtful and also in order to confirm the presence of faint 
lines, spectra obtained with the AURELIE spectrograph 
of the Haute Provence Observatory in July 1997, were 
examined. No variatiations were found, at least for the 
lines studied. The spectral resolution of these AURELIE 
spectra is 10,000 and the signal to noise ratio similar to 
that of the spectrum used in the present work. 

The main requirement of the SAC method to be ap- 
plied is the presence a large number of emission lines of 
different multiplets of the same ion which form in the same 
region of the stellar envelope. P Cygni shows an optical 
spectrum that is quite rich in emission lines (permitted 
and forbidden), including in particular many lines of N 
II and Fc III. The permitted transitions from less excited 
levels are accompanied by blueshifted absorption compo- 
nents, while the lines from more excited levels arc, like the 



forbidden lines, purely in emission (Stahl et al. 1991, Stahl 
et al. 1993, Markova and de Groot (1997). To increase the 
number of lines in our sample we decided to analyse in ad- 
dition to the permitted lines with a pure emission profile 
also those with a P Cygni-type profile. Thus the sample 
of lines, we were able to analyse by the SAC method, con- 
sists of 30 lines from 7 multiplets of N II and 42 lines from 
11 multiplets of Fe III. The multiplets as well as the dis- 
tribution of the lines by multiplets are given in Columns 
1 and 6 of Table 1. 



To construct the SAC of a given multiplet one needs 
to know the lower-level statistical weight, g, the oscillator 
strength, / and the deredenned emission flux, Fa, of 
each line from this multiplet. In the present analysis we 
have used net emission equivalent widths (i.e. equivalent 
widths corrected for photospheric absorption), which we 
converted into line fluxes using a deredenned continuum 
calculated from B = 5.14 and V = 4.72 and deredenned 
with E(B-V)pcy g = 0.63. In deriving the net emission 
equivalent widths (Wa) of the lines from their observed 
profiles, we were confronted by the following problems: 

(i) possible contamination by photospheric absorption; 

(ii) overlapping between wind absorption and emission in 
the case of lines with a P Cygni-type profile; 

(iii) blending effects between pure emission lines which 
are not well-resolved in the spectrum. 

The correction for photospheric absorption was per- 
formed by subtracting a synthethic spectrum, derived, us- 
ing model calculations of Israelian ( 1995 ), from the ob- 
served spectrum. The spectrum, obtained in this way, was 
then measured for the equivalent widths of the lines. It 
must be noted here, that we were not fully convinced that 
such corrections had to be done, nor whether they were 



valid. This is because Markova et al. (2001) showed that 
the wind of P Cygni is probably sufficiently optically thick 
for the star to form a permanent pseudo-photosphere. If 
this is the case, the observed profiles do not need to be 
corrected for photospheric absorption. However, Markova 



and de Groot (1997) found that the emission peaks of the 
P Cygni profiles of Balmcr, He I and N II lines in their 
1990 spectra, were red-shifted with respect to the systemic 
velocity of the star, V sys = -22 km s _1 , a result that can 
be easily explained if one believes that these peaks are dis- 
torted by photospheric absorption. In addition, Israelian 
( 1995 ) argued that the complex structure of [Fe II] forbid- 
den lines in P Cyg's spectrum is due to blending effects 
caused by weak photospheric lines. The results outlined 
above can be reconciled, if one suggests that the wind 
opacity can change significantly, a conj ectur e that seems 
quite reasonable. In fact Markova et al. ( 2001 ) showed that 
the mass loss rate of P Cygni can vary within 26% on a 
timescale of at least 600 days. In any case the photospheric 
correction has a small effect on our results 

The contribution of wind absorption to the P Cygni- 
type profiles was accounted for by means of two methods, 
both of which turned out to give similar results as far as 
SAC slopes are concerned. The first way we proceeded was 
to measure the flux of the red part of each profile, integrat- 
ing the emission flux redward of the emission peak, and 
multiplying the value obtained by 2. This aproach seems 
to be appropriate for our case since emission components 
are expected to be symmetric with respect to their peak 



intensity (Markova & de Groot 1997). The second way we 
proceeded was to fit various Gaussians to the observed 
P Cygni profiles. It turned out that three Gaussions of 
different FWHM (Full Width at Half Maximum) 
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- two centered at V sys for the emission and one for the 
absorption - are usually enough, as a first approximation, 
to obtain a satisfactory fit. The contribution of absorption 
at velocities of more than V sys in this type of fit is small. 

To separate and measure blends between pure emission 
lines we used Gaussian fits to the profiles. This fitting was 
performed taking into account the velocity and FWHM of 
the lines. Blending was also considered taking into account 
the components's location on the self-absorption curve. 
Two examples of blend separations are shown in Fig. 1. 
The numbers given in the upper right hand corner of the 
figure denote the wavelength of the central position, the 
flux and the Gaussian width (in A) of each component. 
Let it be noted that not only the differences in line flux, 
represented by a dashed line, between the fit (dotted line) 
and the observations (solid line), are very small, but also 
that the central velocities and the FWHM of the three 
N II lines (multiplet 36), as determined by the fit, are in 
perfect agreement with each other as would be expected if 
the lines were formed in the same region of the wind. The 
equivalent widths of the single pure emission lines were 
measured, fitting their profiles with a Gaussian. 



citation potentials from 22.54 to 24.59 eV being called 
"multiplet 705" and those with the corresponding excita- 
tion potentials from 23.61 to 25.68 eV being called "multi- 
plet 756". The results of the least squares solutions, where 
at least three lines are present, are given in Table 1. 

The first thing we see is that the vertical dispersions 
of the points around the linear solutions are in most cases 
not very large in the same multiplet. The RMS (root mean 
square) deviation is less than 0.4 for all multiplets except 
multiplets 756 of Fe III and 46 of N II. Secondly, the slopes 
of the graphs, found by linear least squares solutions, are 
still large for more excited multiplets, indicating that their 
lines are optically thick, even though their blue shifted 
absorption is not detectable. This is still true if only the 
more significant solutions are taken, for which the values 
of log (g/A) of some lines in the same multiplet differ by 
more than 0.5. There are even signs of a slope increase for 
the most excited multiplets. The range of log (g/A) values 
hardly changes between the least and the most excited 
multiplets of N II, while in the case of Fe III, the log 
(.9/A) values are however considerably less for the first 3 
multiplets. The slopes and other relevant information are 
listed in Table 1. 



3. Application of the SAC method to P Cyg 



4. Discussion 



At the beginning of this section we would like to note that 
up to now the SAC method has, in all published papers 
except for conference proceedings, been applied to study 
the emission spectrum of the Fe II ion only. In the present 
study, we shall for the first time try to apply this method 
to the spectrum of Fe III and N II. 

In our present work we have firstly tested that the 
points of the graphs of log ^4- against log(g/A) for each 
multiplet lie near a curve. The slopes of the curves were 
then compared for different multiplets. The graphs with 
linear least square solutions for different multiplets are 
plotted in Figs. 2 and 3 

Among these plots are two for different groups of very 
highly excited Fe III lines, those with lower to upper ex- 



The origin of the emission spectrum of P Cygni was partly 
discussed by Wolf and Stahl (1985) who suggested that 
the upper terms of multiplets 115 and 117 of Fe III were 
pumped by two lines of He I. Pumping by two UV tran- 
sitions in He I (AA 522, 537), followed by downward cas- 
cades, was suggested by Markova and de Groot (1997) to 
explain the Fe III emission lines of multiplets 113, 114, 
118 and 119. In addition, the authors noted that the N II 
multiplets with the highest upper term excitation poten- 
tial (~25 eV) might be populated by dielectronic recom- 
bination, but the temperature of the wind seems to be too 



low for this mechanism to work. Stahl et al (1993) how- 
ever suppose that the highest excitation N II multiplets 
are excited by an unknown pumping mechanism. 




Fig. 2. Plots of log £4- against log(gfX) (SAC) for Felll multiplets 4,5,68,113,114,115,117,118,119 and "756' 
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Table 1. Results from the SAC analysis 
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The optical thickness of the more excited Fe III and 
N II lines without blueshifted absorption components, es- 
tablished by us, is quite surprising, while the indications of 
an increase of SAC slope with increasing excitation, if con- 
firmed in future work, would be in fact counter-intuitive. 
At the same log(<?/A) more excited lines are expected to be 
in any case optically thinner. A conspirancy of differential 
non-LTE effects appears unlikely as an explanation, espe- 
cially when the points for a multiplet do not show large 
deviations from a curve. However we should note, that it 
is hard to rigorously test for such an effect, as both the 
upper and lower levels of the observed optical Fe III and 
N II lines could in principle show deviations from LTE. In 
objects showing many Fe II emission lines, the lower levels 
of strong lines are very often in LTE and in that case one 
only needs to test for effects depending on the upper level. 

Another explanation, which could be suggested, will 
not work either, if the increase of curve slope at the same 
log(<7/A) for more excited multiplets is confirmed in future 
work. One might suppose that the more excited lines are 
only emitted from certain parts of a spherically symmetric 
wind. In that case the absence of blueshifted absorption 
components for excited multiplets might be explained by 
a large source function, so the emission in the lines per 
unit surface area, would be larger than that of the photo- 
sphere. In that case extra emission instead of absorption 
would be present on the blue side of the more excited 
lines, resulting in a blueshift of the mean radial velocity, 



as was in fact observed by Markova and de Groot (1997). 
Such a blueshift could be particularly important if part 
of the receding material was behind and so occulted by 
the photosphere, but the symmetry of the lines, noted by 
Markova and de Groot, would suggest that such occulta- 
tion is small. In any case if the wind has spherical symme- 
try, the less excited lines could not be optically thinner, 
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of the wind, covering the photosphere, with only a small 
optically thicker contribution to their emission from the 
clouds. The physical reason for such a situation is however 
not immediately clear. Heating due to shock waves might 
occur in the clouds, producing extra ionizing radiation, fol- 
lowed by more recombinations and cascades plus pumping 
to extremely excited levels Such a situation might also be 
a way of explaining the presence of the most excited N II 
lines. Theoretical work is required to see whether this is a 
viable option. 

We must also finally point out that effects such as an 
increase in self-absorption curve slope for more excited 
multiplcts, would be hard to find in detailed spectral syn- 
theses, based on spherically symmetric wind models. Semi- 
empirical methods can be quite powerful in searches for 
unknown physical processes, present in relatively complex 
situations. 
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Fig. 3. Plots of log against log(gfX) (SAC) for Nil 
multiplets 3,5,19,20,28,36,46 

because those lines would also be emitted in the region of 
emission of the more excited lines, where they should not 
be optically thinner than those excited lines. 

If we exclude the previous explanations, we may need 
to conclude that the wind of P Cyg is anisotropic. 
Optically thick excited lines with no blue shifted absorp- 
tion components would in that case be mainly formed 
in clouds, which either have a large source function or 
which do not occult much of the photosphere. The opti- 
cally thinner lower excitation lines with blucshiftcd ab- 
sorption, could then be mainly formed in large regions 
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